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NON-LTE RADIATION PROCESSES: 
APPLICATION TO THE SOLAR CORONA 



S. COLLIN 

Observatoire de Paris 
F -92195 Meudon, France 



Abstract. These lectures are intended to present a simple but relatively 
complete description of the theory needed to understand the formation of 
lines in non-local thermodynamical equilibrium (NLTE), without appealing 
to any previous knowledge except a few basics of physics and spectroscopy. 
After recalling elementary notions of radiation transfer, the chapter is fo- 
cussed on the computation of the level populations, the source function, the 
ionization state, and finally the line intensity. An application is made to 
forbidden coronal lines which were observed during eclipses since decades. 



1. Introduction 

Astrophysical media, at least for the fraction we can observe, are often very 
dilute compared to those we are used to on Earth. In the interstellar medium 
for instance, the number density is in average one atom per cm 3 , in HII 
regions and planetary nebulae it is about 10 atoms per cm 3 and in the solar 
corona 10 8 atoms per cm 3 , orders of magnitudes less than those one can 
get in laboratory experiments. Therefore the notion of "thermodynamical 
equilibrium" does not hold, and quite unusual phenomena are taking place, 
such as the emission of intense "forbidden lines" never observed on Earth. 
We will in particular show how to compute the intensity of visible forbidden 
coronal lines, and show how these lines, and others, can be used to get 
physical parameters of the solar corona. 
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Figure 1. The specific intensity. 

2. Basics of Transfer 

2.1. PHOTOMETRIC QUANTITIES 

We shall first define a few basic photometric quantities. Note that CGS 
units will generally be used. 

• The specific intensity I u : The energy dE v crossing a surface of area 
dA in the direction of the normal, in a solid angle dfl, during a time dt, in 
a frequency interval dv, is (cf. Fig. 1): 

dE v = I v dv dn dt dA. (1) 

It can be expressed in Watt m -2 ster -1 Hz -1 , but other units are used 
in the visible or in the X-ray range. It can be also defined per interval of 
wavelength, and according to the relation \v = c (speed of light), one gets 
h = Iuc/X 2 . 

The specific intensity can be used both for the source of radiation or for 
the receptor, and it is constant along a path ray in the vacuum. 

• The mean intensity J v : 

Jv = ~k I Iv dn (2) 

• The flux F v : 

It is the power crossing a unit surface per unit frequency interval, in all 
directions (cf. Fig. 2): 

F v = J I v cosO dQ (3) 

F v is generally expressed in Jansky in the radio and far infrared range: 1 J 
= 1(T 23 in CGS, or 10" 26 in MKS (Watts m" 2 Hz- 1 ). 
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Figure 2. The flux. 




Figure 3. The flux received from a star. 

As an application, we can compute the flux from a uniform sphere ra- 
diating isotropically (not a very good approximation for a star, actually), 
cf. Fig. 3): 

F v = [ I v cos 9 dn = I u dcf) sm9 cos 9 d6, (4) 
J Jo Jo 

where 9 C is the angle under which the sphere is seen from the observer. If 
R is the radius of the star, and r its distance, one gets sin# c = R/r, and 

F u = irI u (R/r) 2 . (5) 

Note that the flux at the surface of the sphere is ttI^. 

2.2. TRANSFER EQUATION 

When a light ray does not propagate in the vacuum, the specific intensity is 
not constant: emission adds energy, and absorption removes energy. There 
is also diffusion, in which the global luminous energy is not changed, but 
can be increased in one direction, and decreased in another. We shall not 
take into account the polarization of the radiation, and the possibility of 
non stationary phenomena. 
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2.2.1. Transfer equation in a non diffusive medium 

The monochromatic emissivity r\ v is defined as the power emitted per unit 
solid angle, per unit frequency interval, per unit volume. The monochro- 
matic absorption coefficient Xu is defined as follows: Xvd v is the power 
absorbed per unit solid angle, per unit frequency interval, by a slab of unit 
length normal to the direction of the propagation. Note that here Xv is 
the inverse of a length. Some people use instead the "opacity coefficient" , 
defined per unit mass, and the emissivity per unit mass and not per unit 
volume. The variation of l v on a path length dl in the direction of the light 
ray is therefore: 

dl v = (-XuL + Vu)d£. (6) 

Let us define: 

dr u = - Xvdl and r v = J ~Xud£, (7) 

t v is the optical depth, which decreases towards the observer. 
The transfer equation writes then: 

^ = Iu - S v , (8) 
dr u 

where S v = r] u /xu is called the source function. 

If there are several emission and absorption processes, they have to be all 
taken into account in the source function and in the absorption coefficient. 
For instance, if at a frequency v a line is superposed onto a continuum, 
then S u = S l j ne + S c p on \ and dr v = dr l v ine + dr™ nt . 

2.2.2. Transfer equation in a diffusive medium 

If there is a diffusion process, it must also be taken into account in the 
transfer equation. The diffusion coefficient is defined like the absorption co- 
efficient: a u I u is the power diffused per unit solid angle, per unit frequency 
interval, by a slab of unit length normal to the direction of the propaga- 
tion. Let us assume that diffusion is coherent (i.e. without any change of 
frequency), and isotropic. As a consequence there is a corresponding emis- 
sion coefficient, which is equal to a u J u , and the transfer equation becomes: 

dl u = [{-xv ~ <t v )I v + (o v Jv + Vu)]dt (9) 

If one defines now an extinction coefficient, r*°* = / dr*°* with <ir*°' = 
— {Xv + o- u )d£, one gets for the transfer equation in presence of diffusion: 

% = h - si", do) 
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where S 1 * * = (a u J u + rj u )/(x v + cr u ). Although it is formally similar to 
Eq. 8, it differs in that the intensity appears directly in the second term, 
so it is an integro-differential equation. Actually, the diffusion process is 
a probabilistic one, similar to a random walk, and one can show that the 
distance that a photon will travel before being absorbed is equal to [Xv{Xv + 
^i/)] -1 ^ 2 ) while in the case of pure absorption it is equal to Xu l ■ 
In a purely diffusing medium, the transfer equation writes: 

In the corona, the diffusion process is due to Thomson scattering by free 
electrons. As a first approximation it can be considered as a coherent and 
isotropic process. The diffusion coefficient o v does not depend on frequency, 
and is equal to ffrJV e , where N e is the number of electrons per unit volume, 
and &t is the Thomson cross section, equal to 6.65 10 -25 cm 2 . This process 
is very important for the continuum in the visible range, as it is responsible 
of the emission of the K corona. As a consequence one can consider that 
the corona is a purely diffusive medium for the continuum in the visible 
range. On the other hand we will show that diffusion is negligible in the 
transfer of the coronal lines, on which we will focus later on, so we will not 
consider it in the following sections. 

2.3. APPLICATION TO A PLANE-PARALLEL MEDIUM 

A common utilisation of the transfer equation concerns a stratified medium, 
where all physical quantities are constant on infinite parallel planes. It is 
the usual approximation made for stellar atmospheres (cf. Fig. 4). Often 
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interstellar clouds or diffuse nebulae are considered also as plane parallel 
media. The optical depth is then defined in the direction of the normal, 
and the transfer equation becomes, for a light ray which makes an angle 9 
with the normal to the planes: 

cos# ^ = J„ - S„, (12) 

The formal solution of this equation is: 

I u exp = - / b v exp (13) 

where = cos 0. 

We can apply this solution to different cases, according to the boundary 
conditions. 



2.3.1. for a stellar atmosphere 

In this case the boundary conditions are: 

• no radiation incident on the surface, or I u (fi < 0, t v = 0) = 0, i.e.: 

I v (p < 0,r„) = P ft exp (-fc^) dt, (14) 
•/ o |M| V A 4 / 

• the radiation at = oo remains finite, i.e.: 




and the intensity at the surface is: 




One gets also the flux emerging from the surface: 



F v (fi>0) = 




= 2vr / S v £ 2 (T v )dT v , (17) 
Jo 

where £ n (x) is the order n integro-exponential: £ n (x) = J °° exp(—ux)u~ n du. 

Although Eqs. 16 and 17 are only formal solutions which require to 
know the variation of S v as a function of depth at all frequencies, they can 
be of some help to understand intuitively two observations. 
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First we see that the intensity is approximately equal to the source 
function at t v = [i, and that the internal layers below do not contribute 
to the radiation. It explains the limb darkening effect at the surface of the 
Sun. According to the Eddington-Barbier relation, the source function in 
the photosphere is proportional to r, so it decreases towards the surface (but 
not necessarily in the chromosphere, cf. P. Heinzel's lectures). Thus when 
we observe the limb (/x — > 0), we are seeing layers close to the "surface" 
where S u is small, so I u is also small. When we observe the center of the 
Sun (fi = 1), we are seeing deeper layers where S u is large, and I v is large 
too. So the center of the disk will appear brighter than the limb. 

Second, in the case of stars, we do not observe the intensity, but the flux, 
which is approximately equal to the source function at t„ = 1. Assuming 
again that S v decreases with height in the photosphere, we can understand 
why we see absorption lines in the stellar spectrum. The absorption coef- 
ficient Xv in a line is larger than in the surrounding continuum, so we are 
seeing in a line the layers close to surface, and in the surrounding continuum 
the deeper layers: the lines are then in absorption. The effect is inverted 
if the source function increases with the height, as it is the case for lines 
formed in the chromosphere, or in extended envelopes of stars: the lines are 
then in emission. 

2.3.2. for a homogeneous slab of finite thickness 

One defines an optical thickness increasing towards the observer, dr u = 
Xv d£. Since the medium is assumed homogeneous, the optical thickness of 
the slab is T v = J Xv d£ = Xv H in the direction of the normal, where H is 
the geometrical thickness. 

The boundary conditions are now: 

— incident intensity I v (r v = 0) = I v q, 

— no incident intensity at r v = T v . 

The solution of the transfer equation is thus (for S u = const.): 



If T v < 1, the slab is called "optically thin"; if T v > 1 it is called "optically 
thick". 

There are three interesting cases: 
— 1. a non emissive layer (a cold cloud in front of an intense source): 



I v (T v ,n) = I u0 exp (-— ) + S u 1 - exp (- — 



\ V J Y V H 




(18) 




(19) 
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— 2. an optically thin layer with T„ C 1 (assuming that [i is not <C 1): 
I v (T v ,jj) = I v0 expf-—^ 



= exp M (20) 

V M / A* 

— 3. an optically thick layer with T„ 3> 1: 

I v {T v ,n)~S v . (21) 

For a finite slab with no incident radiation one can show that F v = 
2ttS v [0.5 - S 3 (T V )], and deduce that: 

— 1. in the optically thin case: F v = 2irS u T v = 2irr] u H, or: L u = Airrj^x 
Volume (where L u is the total power emitted by the slab), 

— 2. in the optically thick case: F u = irS u , or L u = irS u x Surface. 

To summarize, in the optically thin case, one "sees" the emissivity, and the 
power is proportional to the volume. In the optically thick case, one "sees" 
the source function, and the power is proportional to the surface. 

Finally note that in a purely diffusing medium, the same equations 
hold, replacing S v by J„, and t v by the diffusion coefficient. So one gets 
in particular for an optically thin medium (i.e. T^* <C 1) with no incident 
radiation on the line of sight: 

I v = . (22) 

This equation applies to the continuum of the solar corona in the visible 
range (the K corona). For a medium which is absorbing and diffusing, and 
optically thin both for diffusion and for absorption, but where the diffusion 
coefficient is negligible compared to the absorption coefficient, the solution 
of the transfer equation becomes: 

/„ = ^ +H «" . (23) 
A* 

This equation applies to the visible lines emitted by the solar corona (cf. 
later) 

3. Local Thermodynamical Equilibrium (LTE) 

3.1. RECALLING THE LAWS OF THERMODYNAMICAL EQUILIBRIUM 

(TE) 

The thermodynamical equilibrium is the stationary state of an ensemble of 
interacting particles and photons which should be achieved in an infinitely 
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thick medium (called a "Black Body") after an infinite time. Photons and 
particles have then the most probable energy distribution, which corre- 
sponds to microreversibility of all processes. For instance, there are as many 
radiative (resp. collisional) excitations from the level A to the level B of an 
atom, as radiative (resp. collisional) deexcitations from the level B to the 
level A per unit time. 



3.1.1. Energy distribution of photons: the Planck law 
It writes: 



I v = B v — 



2hv 3 



exp 



hu\ 
kf) 



i -i 



-1 



(24) 



— h : Planck constant = 6.6262 10 27 erg s, 

— k : Boltzmann constant = 1.3806 10~ 16 erg KT 1 

— T : temperature. 

Caution: it can also be expressed in units of wavelength, and it writes then: 



2hc 2 
A5 



exp 



\XkTj 



(25) 



which has different shape and position of the maximum. 
The integration over v or A gives: 



B 



I 



By dv 



<tT 4 



7T 



(26) 



where a is the Stefan constant = 5.6696 10 5 erg cm 2 s 1 K 4 . 
There are two limiting cases : 

— for hv <C kT, the Rayleigh- Jeans law (used in radio-astronomy): 

2hv 3 kT 2kT 



B v 



hv 



A 2 ' 



for hv » kT, the Wien law (used in the X-ray range): 

2hv 3 



B„ 



■ exp 



hv 
kf 



(27) 



(28) 



3.1.2. Energy distribution of particles in non quantified levels: 
the Maxwell law 

It gives the number of particles per unit volume whose velocity projected 
on an axis z is between v z and v z + dv z : 

( M X 1 / 2 / Mv 2 \ 
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where ./V is the number of particles per unit volume, and M is the mass of 
the particles. This expression is used e.g. to compute the profile of spectral 
lines broadened by thermal Doppler effect (cf. later). 

Since at TE the velocity distribution is isotropic, the number of particles 
per unit volume whose absolute velocity is between v = (v 2 + v 2 + v 2 ) 1 / 2 
and v + dv is equal to: 

Arf M Y /2 ( M(v 2 x + v 2 + v 2 ) \ 
dN = N {^kT 6XP 2kf dv x dv y dv z 



= N 



( M \ 3 / 2 / Mv 2 \ 9 , 



and the number of particles per unit volume with energy between E = 
Mv 2 /2 and E + dE is: 

iN = N {^) m ^{-w)^ EindE - (31) 

This expression is used for instance to compute the rate of collisional exci- 
tations (cf. later). 

Note that the Maxwellian distribution is normalized (/ dN = TV). 

3.1.3. Energy distribution of particles in quantified levels: 
the Boltzmann law 

Consider two levels of energy E m and E n (E rn > E n ), with E m — E n = 
Enm = hv nm , where v nm is the frequency of the corresponding atomic 
transition. The number of particles per unit volume in each level (called 
the "level population") is given by: 

N m 9m ( E nm 



where g m and g n are the "statistical weights" of the levels m and n (they 
are given by spectroscopic tables, as well as the energies of the levels). 

3.1.4. Distribution of particles in different ionization states: 
the Saha law 

It gives the distribution of particles of a given species in a given ionization 

~Nf~ ~ ~W P I kf) ' { } 

— Nl +1 and N} are the numbers of ions on the fundamental level of states 
ionized i + 1 and i times, per unit volume, 
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— N e is number of electrons, per unit volume, 

— gj +1 and g\ are the statistical weights of the fundamental level of states 
i+1 and i, 

— Xi is the- ionization potential of state i (i.e. the energy needed to extract 
an electron from the fundamental level), 

— m is the electron mass (9.1096 10~ 28 g). 
Using Boltzmann law, one gets: 

Nj+iN e 2U i+1 { 2TrmkT \ 3 / 2 f Xi \ 

— Ni and iVj + i are the total numbers of ions in states ionized i and i + 1 
times, per unit volume, 

— Ui and E/j+i are called the "partition functions": 

U i {T)=Y J 9m «p("^) • (35) 

For small T {E\ m /kT >> 1), Ui(T) = g±. For large T, the summation 
must be made on a large number of levels (depending on the density). 

One can show (cf. books on statistical mechanics) that microreversibility 
is equivalent to Maxwell + Boltzmann + Saha equations. 

At TE the optical thickness is infinite and the intensity is isotropic, so: 

B V (T) = I v (t v — > oo) = S v (36) 

which writes: 

r) v = X u B U (T). (37) 

It is the well known Kirchhoff law, which is used to compute the emissivity 
when the absorption coefficient is known (or inversely). 



3.2. FROM TE TO LTE 

At TE, l v is equal to B U (T) and since T is constant, dl u /dr u = and there 
is no transfer of radiation. 

Let us now consider a state where T varies in the medium, but Maxwell, 
Boltzmann and Saha equations are satisfied locally (cf. Fig. 5). In such 
a state microreversibility is achieved. It is called Local Thermodynamical 
Equilibrium (LTE) . We will show later that this state is reached if the 
density and/or the optical thickness are large. The source function is then 
given by the Planck law (but not I u , which varies in the medium and is the 
solution of the transfer equation). 
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Figure 5. Local Thermodynamical Equilibrium. 

4. Source function of a spectral line in non-LTE (NLTE) 

A spectral line corresponds to the radiative transition between two bound 
levels of a given atom. If LTE is not reached, there is no microreversibility 
between the transitions but simply a stationary equilibrium in which the 
number of ALL processes populating a level is equal to the number of 
ALL processes depopulating this level. The main direct and inverse atomic 
processes populating and depopulating a bound level are: 

— collisional excitations and deexcitations 

— radiative excitations and deexcitations 

— collisional ionizations and recombinations 

— radiative ionizations and recombinations. 

Several other processes can be important, such as dielectronic recom- 
binations and auto ionizations, Auger processes, charge exchanges, but we 
will not consider them here. 

4.1. RADIATIVE EXCITATIONS AND DEEXCITATIONS 

Consider again the two bound levels m and n with populations N m and 
N n , energies E m > E n , E nm = E m - E n . 

4.1.1. Einstein probabilities 

The Einstein probabilities A mn , B mn and B nm are defined as follows: 

— the number of radiative spontaneous transitions from level m to level 
n per unit volume per second is N m A mn . Each corresponds to the 
emission of a photon hv = E nm . 

— the number of radiative transitions from level n to level m per unit 
volume per second, corresponding to the absorption of a photon hv = 
Enmi is N n J u B nm . 

— the number of radiative transitions from level m to level n, induced 
by a photon hv = E nm , per unit volume per second and leading to 
the emission of a second photon hv = E nm , is N m J u B mn . It is called 
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"induced emission" (discovered by Einstein). Note that J u is the mean 
intensity at frequency u = E nm /h. 

One can determine the relations between these coefficients in TE. Mi- 
croreversibility between radiative excitation and deexcitation gives: 

N m (A mn + B mn Jy) = N n B nm J v (38) 

Boltzmann equation between levels m and n, and the Planck law J u = B v 
lead to: 



Ju 



exp( — )-l 



(39) 



<B mn J QnBnm 

Identifying the terms depending and not depending on T (as these re- 
lations should be verified for any T), one gets: 

g n B nm = g m B mn , and = (40) 

These relations are equally valid in non-LTE, since A mn , B mn and B nm 
are atomic parameters and do not depend on the physical conditions of 
the medium. The A mn are given in spectroscopic tables (or equivalently 
the oscillator strength f nm , cf. later). Note that a different definition of 
the Einstein probabilities is sometimes used, with the density of radiation 
instead of the intensity in Eq. 38. 

4.1.2. Relations of the Einstein probabilities with the absorption 
coefficient, the emissivity, and the source function 

Spectral lines are broadened by several processes, corresponding to an ab- 
sorption profile ^(^) and an emission profile ^'(i / ), normalized so that 
/ $>{v)dv = 1 and / ^(u)du = 1. The emissivity then writes: 

Vv = ^ N m A mn V(v), (41) 
and the absorption coefficient: 

(42) 

It contains a negative term of induced emission, (N m B mn ) / (N n B nm ) = 
(N m g n ) I (N n g m ) , which can be larger than unity in the radio range (because 
the exponential term in the Boltzmann equation is close to unity). In this 
case the medium becomes an amplifier instead of an absorbant (MASER 
effect). 
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If = $(z/), one gets: 



q i)v f N m A mn \ / N m B mn , (43) 

Xu V N n B nm ) \ N n B 

m 



and the source function is constant along the line profile. 

The absorption and emission profiles are generally the convolution of a 
Gaussian and a Lorentzian function, the first being due to thermal and tur- 
bulent broadening, the latter due to radiative and collisional damping. Such 
a convolution is called "Voigt profile" . The emission and absorption profiles 
are identical in the case of coherent scattering, or if on the contrary there 
is complete redistribution of frequencies in the absorption-reemission pro- 
cess. It is the case for a majority of lines, in particular for the coronal lines 
considered later in this chapter. Since they correspond to forbidden transi- 
tions, radiative damping is negligible, and since they are formed in a dilute 
medium, collisional damping is also negligible. So the emission/absorption 
profile reduces to a gaussian function. Let us rapidly recall the theory. 

If an atom with a transition Hvq = E nm has a velocity projected on the 
line of sight equal to v z , it absorbs or emits photons propagating towards 
the observer, shifted by Doppler effect to v = vq(1 + v z /c). For an assembly 
of atoms with a Maxwellian distribution of velocities, the number of atoms 
able to absorb or emit at a frequency u, in a unit frequency interval, is 
therefore, according to Eq. 29: 

and one deduces for the absorption/emission profile: 



= — =^ exp 



(45) 



where 



(2kT\ l/2 



VW < 46 » 



is called the "Doppler width" of the line. If the gas has also a micro- 
turbulent velocity Vturb, it should be taken into account in the Doppler 
width: 

"0 (2KT , „ 2 



as turbulence corresponds also to a Gaussian distribution of velocities. 
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The absorption coefficient at a frequency v thus writes: 



1 _ Nm Bmn \ (48) 



and the absorption coefficient at the line center is: 

Y -b^ NB - (\ N ™ B ™A (49) 

Instead of the Einstein probability, one often uses the oscillator strength 
f nm which is related to B nm by: 

_ ?re 2 f 4?r 
mc Alt' 

leading to: 

Ar vre 2 1 / iV m B mn \ 
Xv = N n /„ raT r 1 - — — — • 51 

Numerically, Tie 2 /mc = 0.027 in CGS, g n fnm = 1-5 x W s g m A rnn X 2 (A 
in micron, A mn in s _1 ), so / is of the order of unity for permitted lines, 
and <C 1 for forbidden lines. As an example, for hydrogen: 

2 6 1/1 1 V 3 1 1 

J nrn — 7= ~ J 9 ) 5 ' / 

3\/3tt g n \n z m z J m a n 6 
where g is the so-called "Gaunt factor" , of the order of unity. 



4.2. COLLISIONAL EXCITATIONS AND DEEXCITATIONS 

In passing close to an atom, a perturber P, having a kinetic energy £kinj 
can loose kinetic energy in exciting an atom from the bound level n to the 
bound level m. The kinetic energy ££ in of the perturber after the collision is 
equal to E^in ~~ E nm . In the inverse process, the perturber gains the energy 
E nm in deexciting the atom from the level m to the level n. 

One defines the "collisional excitation rate" C nm and the "collisional 
deexcitation rate" C mn as follows: 

— the number of collisional excitations per second per unit volume is 

N P N n C nm 

— the number of collisional deexcitations per second per unit volume is 

N P N m C mn , 

where Np is the number density of the perturbers. 

The perturbers are most often electrons, except in cold interstellar 
medium, where they can be hydrogen atoms or molecules. Let us call a mn (v) 
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(resp.. cr n m(v)) the cross section for a collisional deexcitation (resp.. a colli- 
sional excitation), v being the relative velocity of the interacting particles. 
Generally the mean kinetic energy of different types of particles (in other 
words their temperature) is the same. The velocity of the electrons is thus 
much larger than the velocity of the (more massive) ions. If the perturbers 
are electrons, the relative velocity v reduces therefore to the velocity of the 
electrons. According to the definition of C mn and C nm , one gets: 

C mn (respt.C nm ) =< a mn v > (respt. < a nm v >) s" 1 cm 3 , 

where the average is taken over the velocity distribution of the perturbers, 
i.e. the Maxwellian distribution of the electrons: 

f°° I mv \ 

^kf-J J a ™^^{-2VF e ) Wdv < (53) 

f°° I mv \ 

J Emn ^ 6XP {'^ J W dv ■ (M) 

where T e is the electron temperature (generally equal to the ion tempera- 
ture). One should note that for C nm (T e ) the integration is performed from 
Emn t° 00 j since the perturber must have a kinetic energy larger than E nm 
to be able to excite the transition. The cross sections are determined exper- 
imentally and/or theoretically as functions of v, then integrated over the 
Maxwellian distribution, and C mn (T e ) is tabulated. 

As for the radiative rates, one can determine the relations between C nm 
and C mn in TE: 

— 1. microreversibility between collisional excitation and deexcitation 
gives: 

N n N e C nm = N rn N e C mn , (55) 

— 2. Boltzmann equation between levels m and n gives therefore: 

r , 9m rl ( E nm \ 

Cnm = C mn exp — - (56) 

9n V kT e J 

where T e is the electron temperature and N e the electron number density. 
This relation is also valid in non-LTE. 
Numerically one has: 



8.6310 6 V.„ n , - , 

Cmn = TJo S CHI 3 57 

Te' 9n 



where Q mn is a dimensionless "collision strength" , given in different tables 
and articles. 
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4.3. PHOTOIONIZATIONS AND RADIATIVE RECOMBINATIONS 

A photon with an energy hv larger than the energy required to extract 
an electron from a bound level n, Xn, can ionize an atom in the state n, 
ejecting an electron with a kinetic energy -E^in- This is a "bound- free" pro- 
cess, because the energy of the ejected electron is not quantified, while 
the excitation between two bound levels is a "bound-bound" process. It 
is called photoionization and corresponds to the absorption of the photon 
hv. Photoionization can take place either from outer valence shells, or from 
inner shells, when the ionizing photons have X-ray energies. The inverse 
process, corresponding to the emission of a photon hv, is the radiative re- 
combination. Finally, as for bound-bound processes, there is also stimulated 
recombination induced by an incident photon. 

Note: we do not discuss here dielectronic recombinations and the inverse 
process, autoionizations. Autoionizations are generally unimportant. 

One defines the rate of photoionization per unit volume per second from 
level n as N n P n , and the rate of radiative recombinations onto level n 
(which can include dielectronic recombinations) as N e Ni + \a n . a n is called 
the recombination coefficient. 

The absorption coefficient corresponding to a photoionization from the 
level n can be computed from quantum mechanics, or determined experi- 
mentally. For hydrogenic ions it writes: 



2 A e 2 c 2 R 2 Z A g 1 1 
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where g is a mean "Gaunt factor" , of order unity (for example it is equal 
0.9 for photoionization from level 2). Note that the absorption coefficient 
is proportional to z/ -3 . The dependence in v is generally different for non 



hydrogenic ions. The factor 



hv 



takes into account induced 



l-exp(-- 

recombinations (for the demonstration, see Mihalas, 1978). 

Xun is equal to zero for hv < Xn, which is called the photoionization 
edge. One uses currently also a vn = Xvn/N n , called the photoionization 
cross section. Its value at the edge, a„ n , writes for hydrogenic ions (not 
taking into account induced recombinations): 



*i/0 n 



7.906 10" 18 ^§ cm 
Z A 



(59) 



This coefficient is sometimes extended to non hydrogenic ions by defining 
an "effective" Z, Z + s, s being of the order of unity, but it is better to use 
results from quantum-mechanical computations now available in several 
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databases, such as obtained by the OPACITY project or the CHIANTI 
database. 

Using the Kirchhoff law rj vn = Xun,B u (T) and the laws of TE one can 
compute the emissivity corresponding to the recombination process. For 
hydrogen, it gives: 



j] v = N e N H+ 



C g h 
T 3 /2 n 3 4vr 



exp 
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hv 



kT 



(60) 



where C = 3.26 lO" 6 in CGS. 

We can now write the photoionization rate from a level n 
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which gives, for hydrogen, 
4vr g 



P„ = 2.8 15 10 29 N, 



47r g f° 
h n 5 J Un 



1 — exp — — — 



hv 
kT e 



47T 

— J u dv 
hv 



^ v j -3 
—r dv .cm : 

V 



The recombination rate onto level n is: 

N e N t+1 a n =%T^dv cm- 3 s" 1 , 
h Ju n v 

which gives for hydrogen: 
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(62) 



(63) 



(64) 



where £\ is the order 1 integro-exponential. The total recombination coef- 
ficient is an = a n- As shown by Eq. 64 it depends only weakly on the 
temperature. It is roughly given by an = 2 10 _11 T~ 1 / 2 $ cm -3 s -1 , where 
<I> varies from 1 to 3 for T varying from 10 5 to 10 3 K. These expressions 
are valid for hydrogenic ions with a multiplication by Z 2 . For other ions, 
published tables give recombination coefficients as functions of the temper- 
ature, or absorption coefficients as functions of the photon energy. 



4.4. COLLISIONAL IONIZATIONS AND RECOMBINATIONS 

An electron can also be ejected from an atom in the energy state n, when 
another free electron is passing close to it (or another perturber, but elec- 
trons are generally more important, as for collisional excitations): this is 
called a collisional ionization. After the collision the first electron has lost an 
amount of kinetic energy equal to the kinetic energy of the second electron 
plus the potential energy %„. The conservation of energy writes therefore: 



^kin,! — £kin,l + ^kin,2 + Xn, 



(65) 
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where £kin,i an d ^Lm 1 are the kinetic energies of the perturber before and 
after the collision, and -Ekin,2 is the kinetic energy of the ejected electron. In 
the inverse process, a free electron induces the recombination of a second 
electron on the state n: this process is called collisional recombination, or 
three body recombination (because three particles are interacting). 

The rate of collisional ionizations from a level n, per unit volume and 
time, and the rate of collisional recombinations on level n, are written 
respectively N n N e C nc and N 2 Ni + \C cn . Note that with this definition C nc 
and C cn do not have the same dimension! Like for collisional excitations, 
the ionization rate is equal to < a nc v >, where the average is taken on the 
velocity distribution integrated from the ionization edge, and a nc (y) is a 
cross section either measured or computed. A general expression valid for 
a Maxwellian distribution and for kT e < Xi is: 

C ln « W- 8 oT^ X ^ eV) exp (--gQ (66) 

where o is here the number of optical electrons and Xi(eV) is expressed in eV. 
As usual one finds the relation between the direct and inverse coefficients 
using the TE laws: 

N- 

' . (67) 



C n i — Ci n 



4.5. STATIONARY EQUILIBRIUM OF AN ATOM (STATISTICAL 
EQUILIBRIUM EQUATIONS) 

For each level n, one can write an equation corresponding to the equality 
of all processes depopulating and populating the level n: 

Em^n N m (Vnn j Jv^vdv + A mn + N e C mn (T e )^j 

+N e N t+1 a n (T e ) + N 2 e N t+1 C cn {Te) (68) 



Bum J Ju^udu + A nm + N e C nm {T e ) + P n + N e C nc {T e )^j 



where A mn = for m < n. 

These equations can be simplified using the relations between the in- 
verse processes shown above. This is generally done replacing also the level 
populations by their ratios to the LTE populations, b n , defined as: 

N n = b n N l+1 N e f [—^) exp (^J , (69) 
where Xn is the ionization potential from level n. 
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This set of equations forms a linear system of n equations and n un- 
known (the level populations), when T e , N e , iVj+i and / J v $> v dv are known. 
If the medium is relatively highly ionized, and since it is electrically neu- 
tral and made mainly of hydrogen, the electrons are provided mainly by 
hydrogen, and N e is equal roughly to the number of hydrogen nuclei, so 
N e srf 1.2iV H + w 1.27Vh (the correction of 20% being due to the contribution 
from helium and heavy elements). If the medium is weakly ionized, the elec- 
trons are provided by elements having a low ionization potential, like C + 
and metals, and one cannot estimate a priori the value of N e . We will see 
later how the iVj+i are determined. T e is known through an energy balance 
equilibrium, and / J u & u dv requires the solution of the transfer equation. 
So the problem is very complex, and involves the sophisticated numerical 
methods described in P. Heinzel's lectures. 

To simplify the discussion we will assume that the atom is reduced to 
two bound levels n and m. This picture is correct for some forbidden lines, 
such as the HI 21 cm transition, for the forbidden coronal lines considered 
later, and generally when two levels in the ground configuration can be 
decoupled from the resonant lines and from the continuum, due to a large 
difference in potential energies. A most frequent case is also that of 3 levels 
in the ground configuration, which are decoupled from the other levels and 
from the continuum. For a two level atom the equilibrium equation reduces 
then to: 

N rn fyBmn J J u $ u dv + A mn + N e C mn (T e )^j = 

N n (B nm Jj v $ v dv+N e C nm {T e )). (70) 
Using Eq. (43) for the source function, one gets: 

1 — exp 



E 

~kf: 



x JJu^udv +eB u (T e ) N e Q 

S (line) = with e = — 

1 + e A mn 

(71) 

/ J V Q V dv is a diffusion term, as it does not correspond to any emission or 
absorption of radiation. From this expression we immediately see that if 
the density is high, e dominates and S (line) =B u (T e ) . In other words, the 
transition is at LTE. We can also note that LTE is reached more easily 
when / J v § v dv is small (if the medium is optically thin, for instance). The 
expression of S is actually more complicated in general because Q(v) / 

In the case of a multi-level atom, it is possible to write the source 
function as: 

= JJ^dv + eB„(T e ) + e> 
y ' 1 + e + e" 
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where e' and e" contain the other transitions, and to solve the system by 
iteration. This method was used in the past for a limited number of levels, 
but more sophisticated methods converging much more rapidly are used 
presently (cf. P. Heinzel's lectures). 

4.6. REMARKS CONCERNING THE NOTION OF TEMPERATURES IN A 
NON-LTE MEDIUM 

In LTE a unique temperature is sufficient to describe the medium. This is 
not the case in non-LTE, so it is common to use several temperatures, ac- 
cording to the process considered. The most important temperatures which 
are introduced are: 

— The "electron temperature" T e . It is the temperature given by the mean 
value of the kinetic energy, assuming a Maxwellian distribution; Eq. 
31 gives thus: 

- =4^ < 73 > 

where < v 2 > is the root mean square velocity. This is the temperature 
intervening in the collisional rates C mn and C nm , C cn and C nc . 

— The "excitation temperature" T exc , defined by: 

— eX P ~17r— > ( 74 ) 



which gives: 



N n -B nm J \ N n B nm 



S„ - [^4^) ( 1 - ^#^V = B U (T CXC ). (75) 



(76) 



According to this equation, the absorption coefficient becomes: 

E 



Xu = t~ N n B nm ${v) 
Ait 



1 — exp — 



A negative excitation temperature correspond to a negative absorption 
coefficient, i.e. to a maser or a laser effect. 
— The "radiation temperature" T r , defined by: 

I v = B v (T r ) . (77) 

In the Rayleigh-Jeans regime, it leads to: 

% = Yk U (78) 
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and the solution of the transfer equation for a finite layer becomes: 



T r (0)exp(-^ + T e: 



1 — exp 

A 4 



(79) 



This equation is commonly used in the radio range. 
The "effective temperature" T e g, defined by: 



F = J ' F v dv = itB(T cS ) = oTt% . (80) 

The "color temperature" T c , defined by: 

dJ v dB v {T c ) 



dv du 



(81) 



For instance in a nebula ionized by a hot star, T c is constant in the 
visible and UV range, and equal to the surface temperature of the 
star. Also the color temperature of the solar corona is that of the 
photosphere in the visible range (owing to diffusion). 

5. Application to the solar corona 

5.1. IONIZATION EQUILIBRIUM IN NON-LTE 

To solve the statistical equilibrium equations we need to know the number 
of ions in a given ionization state. If for each level the rate of collisional 
ionizations is equal to the rate of collisional recombinations, Saha's law 
applies. One can show that it occurs for large values of the density (N e > 
10 16 cm' 3 ) which are not reached in the corona and in nebular media. 
These media are therefore in non-LTE for ionization equilibrium, which is 
given by the stationarity equation: 

rate of (photoionizations + collisional ionizations) from all levels 
=rate of (radiative + collisional recombinations) on all levels. 

Actually the problem is simplified by the fact that in these media the 
populations of excited levels are very small, so only ionizations from the 
ground level should be taken into account. 

Other simplifications take place in nebular media and in the corona. 

a. Nebular media: 

HII regions and planetary nebulae are dilute media photoionized by hot 
stars. Their kinetic or electron temperature is relatively low (10 4 K). One 
can show that in these conditions radiative ionizations dominate over col- 
lisional ionizations, and radiative recombinations dominate over collisional 
recombinations in the ionization equilibrium, which reduces to: 
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photoionizations from the ground level 
= radiative recombinations on all levels. 
The photoionization rate depends on the distance, on the temperature, 
and on the luminosity of the hot star. For more details on this subject, one 
can refer to Osterbrock's book (1989). 

b. Solar corona: 

In the corona T e is of the order of a few 10 6 K, while T r is only 5800 K. 
Therefore photoionizations are negligible with respect to collisional ioniza- 
tions. On the other hand the density is relatively low (N e < 10 10 cm -3 ), 
and collisional recombinations are negligible with respect to radiative re- 
combinations. The ionization equilibrium thus writes: 

collisional ionizations from the ground level 
= radiative recombinations on all levels 

or: 

N e Ni C lc {T e ) = N e N l+1 ai {T) (82) 

where ati(T e ) = J2^=i a n(T e ) , and Cj c (T e ) is the rate of ionizations from 
the ground level. 

It can be written approximately: 



N 



i+l 



Ni 



where Pj(T e ) is a tabulated function, of the order of unity, n is the principal 
quantum number, and xh, Xi are the ionization potentials of hydrogen and 
of the ion i. 

This equilibrium is called coronal equilibrium and holds not only in 
the solar corona, but also in the hot interstellar medium, in the hot gas 
of galaxy clusters, in supernovae remnants, etc. Its particularity is that 
Ni/J2N is a function of only T e . As an example we show the ionization 
equilibrium of iron on Fig. 6. We see that for the kinetic temperature of the 
corona, T e ~ 10 6 K, iron should be in the form of highly ionized species, 
Fe +9 to Fe +13 in particular. This figure shows also that the ionic fraction 
is a rapidly varying function of the temperature, implying that if it can 
be obtained from the observations, the corresponding temperature will be 
accurately determined. 



5.2. THE VISIBLE CORONAL LINES 

Since the beginning of the century many lines of unknown origin have been 
observed in the visible during eclipses. They were attributed to a new el- 
ement "the coronium" (remember the "nebulium" of planetary nebulae!). 
Grotrian (1939) understood that the intense red line at 6374A is due to a 
forbidden transition in the ground configuration of FeX (cf . Fig. 7) . 
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Figure 7. Term diagram of FeX showing 4 resonance lines in the EUV range and the 
forbidden visible line at 6374 A. 



Many other lines of highly ionized elements were then identified in the 
visible coronal spectrum, the most intense being identified by Edlen in 1942 
to the forbidden transition in the 3s fundamental configuration of Fe +13 . 
These observations lead to the discovery of the high temperature of the 
corona. These lines are only visible on the limb, during eclipses, or using a 
coronagraph, as the disk is much too bright in the visible range (cf. Fig. 8). 

It is not the same for extreme UV and X-ray lines of highly ionized 
atoms, which are emitted only by the corona, and can therefore be seen in 
emission on the disk itself because there is practically no underlying con- 
tinuum in these bands. When the extreme ultraviolet and the X-ray bands 
were opened to observation with the launch of rockets and satellites, a num- 
ber of permitted coronal lines due to resonant transitions were discovered. 
As an example Fig. 9 displays an image of the Sun obtained with EIT on 
board of SOHO in the FeX/XI lines at 17lA. It shows clearly strong vari- 
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Figure 8. An eclipse seen in the forbidden FeXIV 5303A line 




Figure 9. An image of the Sun obtained with EIT on board of SOHO in the FeX/XI 
lines at 17lA. 

ations of the line intensity, due to variations of the temperature and of the 
density in the corona (cf. below). 

But here we shall focus on forbidden visible lines, because they are more 
simple to handle than resonant lines, and as an example we will determine 
the intensity of the FeXIV 5303A line emitted by the corona. The theory 
of these lines is established since decades, in particular with the pioneering 
work of Pottasch (1963). 

The transition between levels 2 and 1 in the ground configuration of 
Fe +13 producing this line is represented on Fig. 10. One has first to check 
that it can be decoupled from the rest of the atom and from the continuum. 
This requires to show that the rates of radiative and collisional deexcitations 
from all upper levels onto level 2, and the rates of radiative and collisional 
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Figure 10. The transition FeXIV 5303 A. 



recombinations onto level 2, are negligible compared to N\B\2 J J v ^ v dv 
and N\N e C\2 (populating level 2), and that the rates of photoionization 
and collisional ionization from level 2, and the rates of induced radiative 
and collisional excitations to all upper levels from level 2, are negligible 
compared to N2B21 J J v ^ u dv, N2A21 and iV^-ZVeC^i (depopulating level 2). 
For instance collisional and induced radiative deexcitations from excited 
levels are not important because these levels are not highly populated, 
owing to their large energies, and collisional and radiative excitations to 
excited levels are negligible compared to excitations between the levels in 
the ground configuration, owing to the exponential term in the excitation 
rates. Collisional recombinations are negligible because of the low density. 
Actually this is not true for all visible coronal lines, and in several cases 
cascades from highly excited levels, are important in populating the upper 
level of the forbidden transition. 

Let us consider thus a two-level atom in an ionized medium (i.e. per- 
turbers are electrons), and compare the different terms of Eq. 70. The 
atomic data we need are: the Einstein coefficient A21 of the transition, 
equal to 60 s _1 , the statistical weights g\ and g2, respectively 2 and 4, and 
the "collision strength" f^i (cf. Eq. 57 and Fig. 10), of the order of unity 
(it decreases from 3 at T = 10 5 K to 0.3 at T = 10 7 K). 

To compute the rate of radiative excitation, we need J u . We will not 
determine the optical thickness in the continuum at 5303 A, as it would 
require to know the absorption coefficient in the continuum at this wave- 
length, due to free-free process that we have not discussed. But we know 
from observations that the corona is optically thin in the optical band, as 
we see the underlying disk without any absorption. We infer that the whole 
corona is illuminated by the photospheric radiation field, close to a black 
body of temperature To = 5800K. Since the line is narrow (Ai/p <C fo), 
J v is approximately constant over the line profile. The corona being illumi- 
nated only from one side, one gets: 



where i? sun is the radius of the sun and z the height in the corona, and 
one deduces the radiative excitation rate (using the relations between the 



/ J^du = J u = B u (T r ) = f-^V (84) 

J Z \ Z + -Ksun / 
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Einstein coefficients): 



NiB 12 J v = Nx A 2 i 



cxp 



\kT J 



z + R s 



(85) 



For the basis of the corona it is equal to 0.9Aq s _1 cm -3 , and the induced 
deexcitation rate, N 2 ^- J v , is equal to 0.45iV2 s _1 cm~ 3 , therefore negligible 
with respect to the spontaneous deexcitation rate. 

The exponential term in the collisional excitation rate exp ^— jj^^J is 
roughly equal to unity, as £12 <C kT e , so g\C\2 = QiCn- According to 
Eq. 57 one gets thus for the rate of collisional deexcitation N 2 N e C 2 i ~ 
2 x 10~ 9 N 2 N e s" 1 cm" 3 for T e = 10 6 K, and the dependence on T e is 
small. Comparing this with the radiative deexcitation rate N2A21, we see 
that above a critical density N e ~ 10 11 cm -3 , deexcitation would be domi- 
nated by collisional processes. It is a very high density, and one can admit 
that deexcitation is always dominated by radiative spontaneous transitions. 
For the excitation the critical density is much smaller, ~ 10 9 cm~ 3 at the 
basis of the corona, and it decreases with height, as the radiative rate is 
proportional to i? 2 un /(z + -R sun ) 2 . So excitation is sometimes dominated by 
radiative transitions, sometimes by collisional transitions. 

It is important to note here that we have used 3 different temperatures, 
namely the temperature of the underlying blackbody To, the radiation tem- 
perature T r , and the electron temperature T e . All are different from the 
excitation temperature T exc (cf. below). 

As a rough approximation we can consider the corona outside the disk 

as a finite slab of thickness 2R sun ^J 1 — (l — 7^—) ~ 2VHR smi , where H 

is the scale height of the corona (< R S un), with no incident intensity on the 
line of sight. The source function in the line is given by Eq. 71. But we do 
not need this source function to determine the line flux, as we can check that 
the corona is optically thin at the center of the line (and a fortiori in the 
wings). Eqs. 46, 50, 51 and the relations between the Einstein coefficients, 

give: xo = 2.310- 19 Aq. N x can be written: ^^^y ^ ^ ^e, 

where Nh is the number of hydrogen nuclei per unit volume. The iron 
abundance, N ^ , is equal to 4 10 -5 , and the ratio ^ is equal to 0.83 
in the corona (owing to electrons from helium and heavy elements). Let us 
make the conservative assumption that iron is entirely in the form of Fe +13 , 
and that N\ is much larger than the populations of all the other levels, i.e. 
Ni = iV(Fe +13 ). Both assumptions are reasonable in the region emitting 
the FeXIV line, actually. One finds then xo ~ 10~ 23 iV e , or to ~ 10- 23 A4, 
where M e is the number of electrons on the line of sight per unit surface, in 
the emission region. Assuming that the whole corona is emitting the line 
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(which is clearly an overestimation) one gets J\f e =< N e > 2^HR sun , where 
< N e > is the density averaged on the scale height, N e ~ 5 10 s cm -3 , and 
finally one gets for H ~ 5 10 9 cm: To < 10~ 4 <C 1. 

At this stage it is also possible to show that the diffusion coefficient 
neglected in all this chapter is indeed negligible compared to the absorp- 
tion coefficient in the line. Since diffusion is due to Thomson scatter- 
ing, the diffusion coefficient is equal to axN e = 0.66 10 -24 iV e , and we 
immediately see that it is smaller than xo- However, it is not possible 
to neglect the emissivity due to diffusion in the continuum, as it is the 
only emission process for the continuum underlying and surrounding the 
line. According to Eq. 22 for an optically thin purely diffusing medium, 
/*/ = TT J U ~ t t B u (T )/2 ~ 0.3 W~ 24 N e B u (T ) ~ 31Q- 6 B U (T ). This is 
indeed the intensity of the visible continuum observed at the basis of the 
corona. And according to Eq. 23, the intensity in the line should simply 
be added to the continuum intensity, since the absorption coefficient dom- 
inates on the diffusion one. So finally we see that the computation of the 
line intensity can be completely decoupled from the diffusion process. 

The intensity in the line is therefore equal to / rj l/ (z)dz (as the emissivity 
is strongly dependent on the position, one cannot use the homogeneous 
approximation). It is proportional to the emissivity at the line center and 
to the emissivity integrated on the line profile (cf. Eqs. 41 and 45), which 
writes, using Eq. 43: 

Hue = ^N 2 A 21 = ^N^NeCn + B 12 J V ) (86) 

47T 47T 

N 2 is always smaller than Ni, since when excitation is dominated by col- 
lision, N 2 /N\ = N e C\ 2 /A 2 \ < 1, and when it is dominated by radiation 
N 2 /Ni = B 12 J U /A 21 ~ exp(-E 12 /kT ) < 1. (Note that T exc which is 
defined either by N e C\ 2 / A 2 \ = 2 exp(— Ei 2 /kT exc ), or by B\ 2 J v jA 2 \ = 
2 exp(— Ei 2 /kT exc ), is different from the previously used temperatures.) 
Thus Aq ~ iV(Fe +13 ) ~ 4 lO^A^A^Fe+^VA^Fe)], as assumed before. 

Eq. 86 gives thus, in the two extreme cases of N e much larger and much 
smaller than the critical density of 10 9 cm -3 : 

- N e > 10 9 cm~ 3 : j] v ~ %*C 12 iV e 2 [A^(Fe +13 )/iV(Fe)] 

- N e < 10 9 cm- 3 : j] v ~ %-B 12 J v N e [N(Fe +13 )/N(Fe)}. 

The first case corresponds to dense coronal loops, and the second one is 
close to the quiet corona. 

Let us assume that the density is larger than the critical one (actually 
it is often the case for other forbidden coronal lines). The intensity writes 
thus: 

W~4 10-s £ "^ <-"< Fe+1S ' 



m 
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The intensity, and consequently the flux measured at the Earth, is there- 
fore a function of the density and of the temperature in the region where 
iron is in the form of Fe +13 . But iV(Fe +13 )/./V(Fe) depends so strongly on 
the temperature, that one considers it equal to unity in a small range of 
temperature, and zero outside. Moreover C\2(T e ) is not strongly dependent 
on T e . One can thus approximate the intensity by 



where / E{T e )dT e = J N^dz is called the emission measure. 

Therefore, measuring other visible forbidden coronal lines correspond- 
ing to different ions should allow to map the emission measure as a function 
of the temperature. More sophisticated methods involving inversion proce- 
dures can also be used to avoid the rough approximations made above. 

We have focussed here on a forbidden line, but it is clear that the same 
kind of computation can be performed for other types of lines, such as reso- 
nance UV or X-ray lines. The equations giving the population of the upper 
level of the transition would not necessarily have been as simple as in the 
particular case considered here, and consequently the dependence of the 
emissivity on the temperature and on the density would have been more 
complicated, but the whole procedure would have been similar. Then, by 
comparing forbidden lines with resonance lines of the same ion in the ex- 
treme UV and in the X-ray, or more generally pairs of lines of the same ion, 
one can get powerful diagnostics of the temperature and/or of the density 
in the region where this ion is abundant. Besides, the EUV and X-ray lines 
are observed also on the disk, where the forbidden visible lines cannot be 
observed. For instance SilX and SiX lines which are density sensitive have 
been observed with the coronal diagnostic spectrometer (CDS) on board 
SOHO and have allowed to obtain the density and the temperature in ac- 
tive regions on the limb and on the disk as a function of the radius and 
the position angle. One can guess that many other important results will 
be obtained from the data gathered in this mission. 

Finally, one should not forget that the simple study of the line pro- 
files using lines of different elements allows to determine both the kinetic 
temperature and the turbulent velocity, as shown also in this chapter. 
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